AD-A202  634 


tmr F,IF  mp1 


University  of  Colorado  at  Boulder 


i 


r 


DEPARTMENT  OF 
MECHANICAL  ENGINEERING 


DTIC 


‘  88  0  6S. 


a 


Three  Dimensional  Amplification  of  Seismic 
Waves  by  Alluvial  Valleys  of  Arbitrary  Shape: 
Part  I.  Incident  P  and  SV  Waves 

K.  R.  Khair 
S.  K.  Datta 

University  of  Colorado 


CUMER-88-2 


A.  H.  Shah 

University  of  Manitoba 


July,  1988 


DT|« 

t-i-ECTE 


N°V  1  6  1988 


DT  IC 

ELECTE 
NOV  1  6  1988 


« 


ACKNOWLEDGMENT 


The  work  reported  here  was  supported  in  part  by  a 
grant  from  the  National  Science  Foundation  (#  ECE-8518604) 
and  a  grant  from  the  Office  of  Naval  Research  (#  N00014-86- 
K-  0280) .  Partial  support  was  also  received  from  the  Natural 
Science  and  Engineering  Research  Council  of  Canada  (#  A- 
7988) .  One  of  the  authors  (SKD)  gratefully  acknowledges  the 
support  from  the  University  of  Colorado  in  the  form  of  a 
faculty  fellowship  for  the  academic  year  1986-87  during 
which  part  of  this  work  was  done. 


THREE  DIMENSIONAL  AMPLIFICATION  OF 
SEISMIC  WAVES  BY  ALLUVIAL  VALLEYS 
OF  ARBITRARY  SHAPE: 

PART  I.  INCIDENT  P  AND  SV  WAVES 


by 

K.  R.  Khair 
S.  K.  Datta 

Department  of  Mechanical  Engineering  and  CIRES 
University  of  Colorado,  Boulder,  CO  80309-0427 

and 

A.  H.  Shah 

Department  of  Civil  Engineering,  University  of  Manitoba 
Winnipeg,  CANADA  R3T  2N2 


ABSTRACT 


Amplification  of  ground  motion  in  three  dimensions  due  to  arbitrarily  shaped 
cylindrical  alluvial  valleys,  has  been  investigated  in  this  paper.v  A  hybrid  numerical 
method  combining  integral  representation  and  finite  element  discretization  has  been' 
employed.  The  incident  wave  is  either  a  plane  P  or  SV  wave,  propagating  at  an  arbitrary 
angle  to  the  axis  of  the  cylindrical  valley.  Numerical  results  showing  significant  dynamic 
amplification  are  presented  for  various  angles  of  incidence  and  two  different  geometries 
of  the  valleys.  P :  1  -  ■  ^ 
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INTRODUCTION 


One  of  the  basic  problems  in  engineering  seismology  is  the  influence  of  local  geo¬ 
logical  irregularities  on  the  ground  motion  caused  by  propagating  seismic  waves.  Anal¬ 
ysis  of  many  earthquakes  (  e.g.  Mexico  City  earthquake  1985, Sam  Fernando  earthquake 
1971)  show  that  large  amplification  of  ground  motion  occurs  at  topographic  changes  and 
when  the  material  properties  are  softer  near  the  surface  them  that  of  the  surrounding 
medium.  Experimented  studies  by  Kagami  et  ed.(1982,1985),  King  and  Tucker(  1984), and 
Tucker  emd  King(1984)  have  confirmed  this  amplification  phenomenon.  Analytical  solu¬ 
tions  to  this  problem  have  been  obteuned  only  for  simple  geometries  (e.g.  Trifunac  1973, 
Wong  emd  Trifunac  1974a,  Semchez-Sesma  emd  Roeenblueth  1979).  Anedysis  of  more 
complicated  geometries  requires  numerical  techniques  such  as  finite  difference  emd  finite 
elements.  Leirge  geometrical  dimensions  of  geotechnical  problems  limit  the  applicability 
of  these  methods  Usually  to  some  simple  situations  in  two  dimensoins. 

To  avoid  the  difficulties  associated  with  the  discretization  of  the  entire  domain 
the  boundary  integral  equation  method  has  been  found  to  be  effective  (Sanchez-Sesma 
and  Rosenblueth,  1979;  Dravinski  , 1982a, b,c;  Wong,  1982;  Dravinski,  1983;  Dravinski 
and  Mossessian,  1987a)  for  studying  geotechnical  problems.  In  this  method  only  the 
boundary  of  the  scatterer  is  being  discretized.  The  main  advantage  is  in  the  reduction 
of  the  number  of  unknown  variables.  However,  this  method  can  be  usefully  applied  only 
when  the  material  inside  the  boundary  is  isotropic  and  homogeneous  (Kobayashi,  1983). 

In  this  paper,  we  present  a  hybrid  method  which  combines  the  advantages  of  the 
finite  element  technique  and  the  boundary  integral  method.  In  the  past  it  has  been 
applied  to  study  scattering  of  elastic  waves  in  two  dimensions  (Franssens  and  Lagasse 
1984  and  Shah  et  al.,  1982). 


2 


Here,  we  study  the  scattering  of  elastic  waves  due  to  cylindrical  alluvial  valleys 
when  the  plane  of  propagation  has  an  arbitrary  angle  with  the  axis  of  the  scatterer.  Thus 
the  problem  is  three  dimensional  and  to  our  knowledge  this  has  not  been  investigated 
before. 

Numerical  results  are  presented  for  semi-circular  and  semi-elliptical  cylinder- 
ical  alluvial  valleys.  The  incident  plane  P  and  SV  waves  propagating  at  arbitrary 
angles  are  considered.  Results  for  the  particular  case  when  the  direction  of  incidence 
lies  in  a  plane  perpendicular  to  the  axis  of  the  cylindrical  valley  have  been  found  to 
agree  with  previously  published  results.  It  Is  shown  that  amplification  can  be  larger 
for  directions  of  propagation  oblique  to  the  axis  than  when  they  are  perpendicular. 

STATEMENT  OF  THE  PROBLEM 

Figure  1  shows  the  geometry  of  the  problem.  As  shown  cylindrical  inclusion  of  arbitrary 
surface  shape  with  axis  parallel  to  the  y-axis  (  not  shown  )  is  perfectly  bonded  to  ho¬ 
mogeneous,  isotropic  elastic  half-space.  Let  us  consider  two  artificial  boundaries  C  and 
B.  The  medium  is,  now,  divided  into  two  regions.  The  interior  region  Ri  is  bounded  by 
the  boundary  B  and  part  of  the  free  surface.  The  exterior  region  Ra  is  bounded  by  the 
free  surface  and  the  contour  C  and  extends  to  infinity  in  the  x  and  z  directions.  The 
area  between  the  contours  C  and  B  is  shared  by  both  regions.  All  the  inhomogeneities 
are  assumed  to  lie  wthin  C. 

We  consider  plane  harmonic  P  or  SV  wave  incident  at  an  arbitrary  angle  to  the 
y-axis  (  see  Fig.  2).  Thus  the  displacement  has  all  three  components. 

Let  Uj  be  the  displacement  component  in  the  i-th  direction  in  a  Cartesian  frame 
and  let  Ta  be  the  second  order  Cauchy  stress  tensor  (ij  =  1,2,3)  having  time  harmonic 
behavior  of  the  form  e-,w*  where,  uj,  is  the  circular  frequency  (rad/sec).  In  each  region 
Uj  and  Ta  satisfy  the  equation  of  motion  given  by  Eq.  (1), 

Tij.j  +  pv2“i  =  ~fu  (t,;  =  1,2,3)  (1) 
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where  p  is  the  mass  density,  fi  is  the  force  per  unit  volume,  and  the  factor,  e~,wt,  has 
been  dropped.  The  stress  free  boundary  conditions  along  the  surface  of  the  half-space 
are  specified  by, 


O  zz  =  0 


erz,  =  0 


(T  zy  —  0. 


The  scattered  field  also  satisfies  the  radiation  condition  at  infinity. 

(I)  The  Exterior  Region  R*, 

In  this  region  the  displacement  Uj  is  composed  of  two  parts, 

ii^tt^  +  u^,  (2) 

where  is  the  free-field  displacement  (including  the  incident  waves  and  their  reflec¬ 
tions  from  the  free  surface)  and  is  the  contribution  of  the  scattered  waves.  The 
scattered  displacement  field  is  represented  by  a  surface  integral  as  discussed  below. 

Starting  with  Betti’s  reciprocity  theorem  (  Aki.fr  Richards  1980),  a  pair  of  solu¬ 
tions  to  Eq.(l)  for  the  displacements  can  be  written  as 


J  J(L’  £_£'2i)  =  y 


where  use  has  been  made  of  the  radiation  condition.  In  the  above  u,  t  represent  the 
displacement  and  surface  traction  caused  by  body  forces  /,  while  v,  s  are  the  displace¬ 
ment  and  the  surface  traction  due  to  body  forces  g  in  region  R0.  The  scattered  field  is 
taken  to  be  the  first  field.  The  second  field  is  the  Green’s  solution.  The  scattered  field 
has  no  sources  in  R0,  hence  /  =  0  .  For  the  Green’s  displacement  field  the  source  is 


represented  by 


9  =  *(L-r >-<"<+**’»ei, 


where  |(r  —  r')|  =  \/(z  —  x1)2  +  (z  —  z')2  ,  e^  is  the  unit  vector  in  the  i-th  direction,  and 
ki  =  uz/cj,  ci  being  the  shear  wave  speed  in  R0.  This  represents  a  line  source  at  (x',  z') 


varying  in  the  y-direction  with  wavelength  The  expression  for  the  Green’s  dis¬ 
placement  tensor,  Gtj,  and  the  corresponding  stresses,  are  presented  in  Appendix 
A.  Substituting  the  expression  for  g  into  Eq.(3)  and  using  the  divergence  theorm  , yields, 
dropping  the  factor  e'*kiy, 

"**•  (4) 

c 

In  writing  Eq.(4)  the  condition  of  vanishing  traction  on  the  free  surface  z=0  and 
the  radiation  condition  at  infinity  satisfied  by  the  scattered  and  the  Green’s  displacement 
fields  have  been  utilized.  The  contour  integration  is  carried  out  in  clockwise  direction 
and  nk  sure  the  components  of  the  outward  unit  normal  vector  to  C.  Here  we  choose 
the  source  points  ( x ',  z')  to  be  on  the  contour  B  inside  Ro  (Figure  1).  This  precludes 
the  singularities  arising  from  evaluating  the  integral  with  the  source  and  the  observer 
at  the  same  location. 

Now  we  apply  Eq.  (3)  to  the  part  of  region  R\,  which  includes  the  scatterer  and 
bounded  by  contour  C.  The  two  fields  are  the  Green’s  field,  with  its  source  outside  this 
region,  and  the  incident  field.  The  incident  field  has  no  sources  here.  This  yields 

0  =  fiGitfi  -ufVijk)  (—nit)  dc.  (5) 

c 

The  integration  along  C  is  calculated  in  counterclockwise  direction,  and 
denote  the  displacement  and  stress  tensor  associated  with  the  incident  field,  respectively. 
Reversing  the  direction  of  integration  of  Eq.  (5)  and  adding  it  to  Eq.  (4),  using  Eq. 
(2),  yields 

Ui(x\  z')  =  £ (GijTjk  -  UjZijk)nkdc  +  uJo),  (6) 

c 

Eq.  (6)  is  the  integral  representation  of  the  total  field  at  any  point  in  the  exterior  region 

R0. 
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(II)  The  Interior  Region  Ri 


This  region  encloses  all  the  inhomogeneities  and  anisotropies.  In  order  to  get 
the  solution  in  region  R\  ,  we  use  the  finite  element  technique.  In  this  approach,  the 
area  of  interest  R j  is  divided  into  a  number  of  elements.  In  each  element  the  particle 
displacement  is  given  by 

«.  =  {*/}T{u«}  (7) 

where 

i  =  1 , 2,  3 

l  =  1, . n  n  is  the  number  of  nodes  per  element, 

and  {$;}  is  an  nxl  column  vector  representing  the  shape  functions  for  the  ele¬ 
ment, 


{*,}  =  [#!  (8) 

In  Eq.(7)  {uj;/}  is  the  particle  displacement  vector  corresponding  to  the  nodal  points 
of  each  element,  T,  {•},[•]  denote  a  transpose,  a  column  vector,  and  a  row  vector, 
respectively.  The  strain  within  an  element  related  to  the  displacement  field  {u}e  is 
given  by 


{<}.  =  [£][*]{«}., 
=  (B]{u}„ 

where  [  L  ]  is  an  operator  matrix 


r  9 
37 

0 

0  ' 

0 

if 

0 

(£]  = 

0 
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0 
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and 


[B|  =  [I]|*l  (12) 

In  order  to  determine  the  elemental  impedance  matrix,  let  us  consider  the  energy  func¬ 
tional  E, 


E  =  1/2  ///<-  .e*  —  pu2  u.u*)dxdydz  —  1/2  J  J(tB.umB  +  tg.us)ds,  (13) 
where  T  and  e  are  stresse  and  strain  column  vectors,  respectively,  defined  as 


T  =  {T}  =  [TxxT„TxxTtyTxxTxv}T 


and 


«  =  {«}  = 

represent,  respectively,  the  traction  force  and  the  displacement  at  the 
boundary  B.  For  the  purpose  of  eliminating  the  intergration  over  the  y-direction,  it  can 
be  shown  that  if  one  takes  the  integration  over  one  wavelength  A,  Eq.(13)  becomes 

E/A  =  1/2  J  JiT.e*  -  pu)2u.um)dxdz  -  1/2  J +  t^.UB)dc  (14). 

The  stresses  T  are  related  to  the  strains  e  by  the  constitutive  equation 

{T}  =  [D]{e}.  (15) 


Using  Eq.(15)  and  Eq.(7)  in  Eq.(14)  and  taking  the  variation,  the  equation  of  motion 
for  region  Ri  can  be  written  as 


$//  Sib  1  /  \  /  *7  1 

.  Sbi  SbbJ\ubJ  \  Fa  J 


(16) 
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where  the  elemental  impedence  matrices  S/j  are  represented  by 

[S|«  =  /  /([B-iriDlffll,  -  pwJ[$,|T[*.|)<k<fe  (17) 

Ac 

and  the  nodal  force  vector  due  to  applied  surface  traction, 

{Y),=  J[B-)Jm*,}{n}dc.  (18) 

c. 

Here  Ae  and  Ce  denotes,  respectively,  the  area  and  the  boundary  of  an  element  and 
{n}  is  the  unit  vector  normal  to  Ct.  [  D  ]  is  6  x  6  matrix  of  the  element  material 
elastic  constants.  It  is  clear  from  Eq.(17)  that  [S]e  is  hermitian  matrix.  In  Eq.(16) 
{ui}  and  {ug}  represent  the  interior  and  boundary  nodal  displacements,  respectively. 
{Yb}  represent  the  interaction  forces  between  regions  Ro  and  f?i  at  the  boundary  nodes. 
Since  there  are  no  forces  on  the  interior  nodes,  hence  {Yj}  =  0,  and  Eq.(16)  becomes, 

s'b'i  Sbb] {us}  =  {y„}  (l9) 

Thus  by  using  the  upper  part  of  Eq.(19),  the  boundary  nodal  displacement  {us} 
can  be  related  to  the  interior  nodal  displacement  {u/}  as 

{u/}  =  -[5//]-1[5/s]{ub}-  (20) 

Now  ,  Eq.(6)  and  Eq.(19)  should  be  combined  in  order  to  solve  for  the  displacement  at 
the  boundary  nodes  Nb-  To  achieve  this  Eq.(6)  should  be  evaluated  at  the  nodes  on  the 
boundary  B. 

Using  the  constitutive  relation  to  express  Tij  and  Eq.(7)  to  present  the  displace¬ 
ment,  Eq.(6)  becomes 

■<*',*')  ««t*>(*\*')+  fmmBr]  -  [*;lr[E)){n)Jc  {u,J 

.c 

+  j{[G\[D\\BB]  -  [*BlT|E))(n)rfc  {uB),  (21) 

.c 
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where  [ Bi ]  =  [£][$/]  and  [Bb]  =  [£][$b].  Evaluating  the  integrals  at  all  the 
nodes  Nb  on  the  boundary  B,  Eq.(21)  can  be  written  as 

{us}  =  [i4B/]{u/}  +  [Abb]{ub}  +  {u^},  (22) 

wherr  [As/]  is  3 Nb  x  3 Ni  and  [Abb]  is  3 Nb  x  3 Nb  complex  matrices.  Substituting 
Eq.  (20)  into  Eq.  (22),  and  collecting  cerms,  yields 

{ub}  =  -{^b/](5//]"1[S/b]{ub}  +  [Abb]{«b}  +  {ub*}  (23) 

Once  {us}  is  solved  from  Eq.  (23)  we  then  vise  Eq.  (20)  to  find  {t»/}.  In  the  following 
we  present  numerical  results  for  the  displacement  amplitudes  on  the  surface  of  the  valley. 

Numerical  Results 

Normalization  of  Variables. 

All  the  variables  used  in  this  paper  are  presented  in  dimensionless  form.  The 
material  properties  of  the  valley  are  normalized  with  respect  to  those  of  the  half  space. 
For  that  purpose  we  have  chosen  the  shear  velocity  /30  and  the  shear  modulus  /i0  to  be 
unity,  distances  are  normalized  with  respect  to  half  the  width  of  the  finite  element  region 
H,  and  the  displacements  are  normalized  to  the  absolute  amplitude  of  the  displacement 
of  the  free-field  motion  ( y| Ui°^l2  +  |u^0  |2  +  |u«°*|2).  The  normalized  frequency  <  is 
defined  as 

k2H 

where  k2  is  the  shear  wave  number  in  the  half  space.  The  poisson’s  ratio  for  all  the 
materials  is  taken  to  be  1/3. 

Testing  the  method 

Sufficiently  general  computer  code  has  been  written  to  investigate  the  three  di¬ 
mensional  problem.  In  order  to  gauge  the  accuracy  of  the  results  we,  first,  let  the 


9 


material  of  the  alluvial  valley  to  have  the  same  properties  of  the  Half  space.  The  cal¬ 
culated  displacements  should  be  equal  to  the  free  field  displacements.  Deviation  from 
the  free  field  response  is  taken  as  a  test  of  the  inaccuracy.  When  the  error  was  found  to 
be  unacceptable  we  increased  the  number  of  elements  in  the  interior  region  R\ .  It  was 
found  that  6  elements  per  wavelength  are  suffecient  to  achieve  the  required  accuracy. 

Second,  the  generality  of  the  program  was  tested  by  running  the  problem  of 
scattering  of  plane  harmonic  P  and  SV  waves  by  semi-circular  alluvial  valley  when 
£  =  0.  This  is  a  plane  strain  problem  and  was  studied  by  Dravinski  and  Moesessian 
(1987a)  for  weakly  inelastic  materials.  Figures  3  and  4  show  very  good  agreement  of 
the  results  for  all  the  angles  of  incidence  and  both  types  of  incident  waves  (  P  and  SV 
)  and  e  =  n. 

Semi  —  circular  alluvial  valley. 

For  semi-circular  valley  the  finite  element  part  is  compsed  of  190  elements. 
Isoparametric  elements  of  3  and  4  nodes  are  used.  The  number  of  boundary  nodes 
Ng  is  29.  The  number  of  the  interior  nodes  IV/  is  178.  The  depth  of  the  finite  element 
region  in  the  z-direction  is  taken  to  be  equal  to  H.  The  results  presented  here  are  for 
e  =  7r  .  The  area  between  contours  B  and  C  is  occupied  by  one  layer  of  elements.  The 
shear  modulus  of  the  valley  fii  is  taken  to  be  1/6  and  the  shear  velocity  0i  is  1/2.  Fig¬ 
ures  5  and  7  show  the  normalized  amplitudes  of  ux,ut,  and  ux  when  P  and  SV  waves 
are  incident  in  a  plane  making  an  angle  45  degree  with  the  axis  of  the  valley.  Figures  6 
and  8  show  the  corresponding  results  when  the  plane  of  propagation  is  parallel  to  the 
yz  plane.  Figure  5  shows  large  amplification  in  the  y-component  of  the  ground  motion 
when  P  wave  propagates  with  9\  =  60°.  Also,  figure  7  shows  large  amplification  when 
the  angle  of  incidence  of  SV  wave,  02 ,  is  30  degree.  From  figures  5  and  7  it  is  seen  that 
uy  is  generally  amplified  the  most  compared  with  ux  and  uz.  Comparison  of  Fig.  5-8 
shows  that,  exept  when  SV  wave  is  incident  in  the  plane  of  yz,  u9  is  generally  amplified 
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more.  For  SV  wave  incident  in  the  yz-plane  one  finds  large  amplification  in  vertical 
displacement  when  02  =  60°.  Also  of  interest  is  to  note  that  larger  amplification  occur 
for  incidence  in  the  45°  plane. 

Semi  —  elliptical  alluvial  valley 

In  this  case  different  finite  element  mesh  is  used.  The  total  number  of  elements 
is  146  ,  Nb  and  iV/  are  equal  to  29  and  131,  respectively.  The  ratio  of  the  minor  axis 
to  the  major  axis  is  0.7  and  the  major  axis  is  along  the  x-axis.  The  material  of  the 
semi-elliptical  valley  has  the  same  properties  as  that  used  in  the  semi-circular  case,  e  is 
again  taken  to  be  tc.  Figures  9  and  10  demonstrate  the  results  for  the  plane  strain  case 
while  figures  11,  12,  13,  and  14  show  the  results  for  different  angles  of  incidence.  It  is 
seen  that  larger  amplifications  occur  in  a  semi-circular  valley  than  in  a  semi-elliptical 
one.  Also,  it  is  found  that  for  P-waves  amplification  increase  with  0\ ,  but  for  S V-waves 
they  remain  about  the  same  first  as  increases  and  then  decrease.  Here  the  most 
amplification  of  displacements  occurs  when  P  wave  propagates  in  a  plane  parallel  to  the 

axis  of  valley  and  $i  =  60* and  when  S V  wave  propagates  in  that  plane  but  with  angle 

•a 

02  =  30.  In  contrast  with  the  semi-circular  valley  we  find  that  uy  is  amplified  most  for 
the  semi-elliptical  valley. 

Conclusion 

A  hybrid  numerical  method  hes  been  employed  to  investigate  the  three  dimen¬ 
sional  amplification  of  ground  motion  due  to  arbitrarily  shaped  cylindrical  soft  valley. 
This  method  combines  the  integral  representation  and  the  finite  element  techniques. 
Numerical  results  are  presented  for  different  angles  of  incidence.  These  results  show 
large  local  amplification  when  the  plane  of  propagation  is  oblique  to  or  contains  the 
axis  of  the  valley.  This  finding  is  rather  important,  because  so  far  only  two-dimensional 
problems  have  been  studied.  It  also  shows  significant  influence  of  the  shapes  of  the 
valley  on  the  ground  motion. 
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This  method  is  also  applicable  to  study  the  scattering  of  SH  and  Rayleigh  waves 
and  to  a  layered  medium.  These  results  will  be  communicated  later.  The  very  important 
feature  of  this  method  is  that  the  calculations  of  Green’s  function  and  stresses  are 
independent  of  the  shape  and  the  material  of  the  valley.  This  means  that  for  different 
shapes  of  valley  and  the  same  frequency  the  Green’s  functions  and  stresses  has  to  be 
calculated  only  one  time.  This  is  the  most  time  consuming  part  of  the  solution.  This 
differs  from  an  indirect  boundary  integral  equation  formulated  on  the  boundary  of  the 
scatterer  for  which  Green’s  tensor  has  to  be  recalculated  at  each  change  in  the  geometry 
of  the  scatterer. 
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APPENDIX  A 


Green’s  Function  and  Stresses. 


(I)  Displacements 
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ki  —  Longitudinal  wave  number  =  u/ci . 

Ci  =  Longitudinal  wave  velocity  in  Ra- 

Note  that  in  the  above  the  coordinates  x,y,z  have  been  nondimensionalized  with 
respect  to  H  and  the  ±  sign  correspond  to  z  >  or  <  z\  respectively. 
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APPENDIX  B 


Incident  Displacement 


(I)  Pwave 


Let  us  consider  P-wave  incident  in  the  XZ  plane.  The  incident  and  reflected 
potentials  can  be  written  as  : 
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The  components  of  the  displacement  vector  in  the  XZ  plane,  Ux,Uz »  can  be 
derived  from. 
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where  $  =  and  ¥  =  . 
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Displacements  in  the  xz  plane  are  related  to  those  in  the  XZ  plane  by  the  following 


relations, 


ux  =  Uxcos<t>  and  it,  =  Uxsin< <f>. 


The  full  expressions  for  ux,u9,  u*  are  documented  below. 
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<*i  =  «7i~ ,  <*2=  *72*,  a3  = -e ^2 


k  =  T]  sind  1  co3<j> ,  £  =  tj  sinffi  sin<f> ,  A'2  =  £2  +  fc2 


=  angle  of  incidence  with  the  negative  Z-axis  in  the  XZ  plane  (0°  <  9X  <  90°). 


<f>  =  angle  between  xz  plane  and  XZ  plane  (0®  <  ^  <  180®  ). 


(II)  SVwave 


Let  us  consider  incident  SV-wave  in  the  XZ  plane: 
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Following  the  same  steps  as  before  in  the  P-wave,  we  get: 
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and  &2  is  the  incidence  angle  with  the  negative  Z-axis  in  the  XZ  plane. 
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List  of  Figures 

Fig.l  :  Geomtry  of  the  problem. 

Fig.2  :  Incidence  angle  in  3-D  defined  in  terms  of  <f>  and  9. 

Fig.3  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Circular  Alluvial  Valley. 
(e  =  7r ,  fj.0  =  l,  v0  =  1/3  ,/?„  =  1.0 
Hi  =  1/6,  i/i  =  1/3, 0i  =  0.5  ) 

(a)  Incident  P-wave:  (f>  =  0°  ,  9X  =  0°. 

(b)  Incident  P-wave:  <f>  =  0°  ,  9i  =  30°. 

(c)  Incident  P-wave:  <j>  =  0°  ,  9\  =  60°. 

Fig.4  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Circular  Alluvial  Valley. 

(a)  Incident  SV-wave:  <f>  =  0°  ,  02  =  0°. 

(b)  Incident  SV-wave:  <j>  =  0°  ,  d2  —  30°. 

(c)  Incident  SV-wave:  <f>  =  0°  ,  92  =  60°. 

Fig.5  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Circular  Alluvial  Valley. 

(a)  Incident  P-wave:  <f>  =  45°  ,  9\  =  30°. 

(b)  Incident  P-wave:  <f>  =  45°  ,  9\  —  60°. 

Fig.6  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Circular  Alluvial  Valley. 

(a)  Incident  P-wave:  <f>  =  90°  ,  0\  =  30°. 

(b)  Incident  P-wave:  <f>  =  90°  ,  9y  =  60°. 

Fig.7  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Circular  Alluvial  Valley. 

(a)  Incident  SV-wave:  <f>  =  45°  ,  02  =  30°. 

(b)  Incident  SV-wave:  4>  =  45°  ,  02  =  60". 

Fig.8  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Circular  Alluvial  Valley. 

(a)  Incident  SV-wave:  <f>  =  90°  ,  02  =  30". 

(b)  Incident  SV-wave:  <f>  =  90"  ,  92  =  60". 
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Fig.9  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Elliptical  Alluvial  Valley. 
(e=  7T,^0  =  1,  v0  =  1/3, /?0  =  1.0 

Hi  =  1/6,  vx  =  1/3, 0i  =  0.5,  Minor  Axis/Major  Axis  =  0.7  ) 

(a)  Incident  P-wave:  <f>  =  0°  ,  6X  =  0°. 

(b)  Incident  P-wave:  <f>  =  0°  ,  =  30°. 

(c)  Incident  P-wave:  <(>  =  0°  ,  =  60°. 

Fig.  10  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Elliptical  Alluvial  Valley. 

(a)  Incident  SV-wave:  <f>  =  0°  ,  02  =  0°. 

(b)  Incident  SV-wave:  <f>  =  0°  ,  02  =  30°. 

(c)  Incident  SV-wave:  <f>  =  0°  ,  02  =  60°. 

Fig.  11 :  Normalized  Surface  Displacement  Amplitude  for  Semi-Elliptical  Alluvial  Valley. 

(a)  Incident  P-wave:  <f>  =  45°  ,  Q\  =  30°. 

(b)  Incident  P-wave:  <f>  =  45°  ,  9X  '=  60°. 

Fig.12  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Elliptical  Alluvial  Valley. 

(a)  Incident  P-wave:  <f>  =  90°  ,  6X  =  30°. 

(b)  Incident  P-wave:  <f>  =  90®  ,  6X  =  60®. 

Fig.  13  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Elliptical  Alluvial  Valley. 

(a)  Incident  SV-wave:  (f>  =  45®  ,  02  =  30®. 

(b)  Incident  SV-wave:  <f>  =  45®  ,  02  =  60®. 

Fig.14  :  Normalized  Surface  Displacement  Amplitude  for  Semi-Elliptical  Alluvial  Valley. 

(a)  Incident  SV-wave:  <f>  =  90®  ,  02  =  30®. 

(b)  Incident  SV-wave:  4>  =  90®  ,  02  =  60®. 


2 


z 

Fig.  1  Geometry  of  the  problem 


Fig.  2  The  arrow  indicates  the  incident  plane  wave 
direction  which  lies  in  the  plane  of  Xz  that 
makes  an  angle  (f>  with  the  plane  of  xz 
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